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Abstract 

Using Fermat's principle, we analyze the effects of very long wavelength 
gravitational waves upon the images of a gravitationally lensed quasar. 
We show that the lens equation in the presence of gravity waves is equiv- 
alent to that of a lens with different alignment between source, deflector, 
and observer in the absence of gravity waves. Contrary to a recent claim, 
we conclude that measurements of time delays in gravitational lenses 
cannot serve as a method to detect or constrain a stochastic background 
of gravitational waves of cosmological wavelengths, because the wave- 
induced time delay is observationally indistinguishable from an intrinsic 
time delay due to the lens geometry. 
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I. INTRODUCTION 



A stochastic background of gravitational waves of cosmological wavelengths may arise 
in the early Universe, for instance as a consequence of quantum effects during a period 
of inflationary expansion, or as the result of gravitational radiation by oscillating cosmic 
strings. Its presence could be manifested as a large angular scale anisotropy in the cosmic 
microwave background, induced by the Sachs- Wolfe effect the differential redshifting 
of photons in the presence of tensor metric perturbations. It is possible that a significant 
fraction of the anisotropy measured by the COBE DMR experiment [[J is due to cosmological 
gravitational waves. [§] So far, the microwave observations cannot determine how much of 
the anisotropy is due to tensor perturbations (gravitational waves) and how much to scalar 
(energy-density) fluctuations. 

Another potential method to reveal the presence of gravitational waves of cosmological 
wavelengths was recently suggested by Allen, |||| namely, to use measured time delays 
between gravitationally lensed multiple images of distant quasars. Gravitationally lensed 
multiple images of a source such as a quasar arrive at the Earth at different times if the 
source, deflector (the lensing body), and observer are not in perfect alignment, because 
there is a difference in geometric path lengths between, and in the deflector's gravitational 
potential traversed by, the different light rays. We shall call these two effects the 'intrinsic' 
time delay of the lens. For a lens geometry where L is the distance between observer 
and deflector and 2r] is the angular separation between the images, the typical intrinsic 
time delay is AT ~ Lrf. Any 'extrinsic' perturbation to the spacetime metric (i.e., not 
associated with the lens itself) would be expected to cause additional time delays between 
the images. For example, an additional time delay would be induced by a cosmological 
background of very long wavelength gravitational waves. 0,|5| As shown by Allen, for waves 
with frequency uj ~ L , the gravity wave-induced time delay is of order AT Lhrj, where 
h is the dimensionless amplitude of the gravitational wave. Therefore, waves of amplitude 
h > r] would be expected to have drastic effects on lens time delays. 
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Based on this effect, Allen claimed that gravitational lenses could serve as gravitational 
wave detectors, j|,|5| and that a bound could be placed on the amplitude of gravitational 
waves of cosmological wavelengths from the requirement that the wave-induced time delay 
in the double quasar 0957+561 not exceed the observed delay of 1.48 yr. |6| (The observed 
delay is generally attributed to the intrinsic delay of the lens.) For 0957+561, the image 
angular separation is 3 x 10~ 5 rad = 6.1 arcsec, and Allen obtained the bound h < 2 x 10~ 5 . 
A gravity-wave background which nearly saturates the above bound would have important 
implications for gravitational lens models and would seriously compromise attempts to use 
lens time delays to measure the Hubble parameter. 

Subsequent to Allen's work, the microwave anisotropy bound on the amplitude of cosmo- 
logical gravitational waves has been significantly tightened by COBE. Through the Sachs- 
Wolfe effect, gravitational waves induce a temperature anisotropy of order their dimen- 
sionless amplitude. From the COBE detection of the quadrupole anisotropy, it follows 
that h ^$ (5T/T)£ = 2 ~ 6 x 10 -6 for wavelengths comparable to the present Hubble radius, 
A ~ Hq 1 = 3000 h" 1 Mpc. This bound is roughly a factor of three smaller than Allen's 
limit. However, although 0957+561 is the first gravitational lens system for which a time 
delay has been reliably measured, other lens systems are also being monitored; in particular, 
for a lens with smaller image angular separation rj and thus smaller intrinsic time delay, the 
wave-induced delay would be even more important, and the corresponding lens bound on 
h potentially more restrictive. Inflationary models suggest that a significant fraction of the 
quadrupole anisotropy could be due to gravitational waves [|],[7| . If this is the case, then time 
delays induced by gravitational waves in gravitationally lensed quasars would be significant. 

In this paper, we reconsider Allen's proposal. Our central theme is that, for measure- 
ments of time delays in gravitational lenses to serve as gravitational wave detectors, the 
observer must be able to separate the wave-induced time delay from the intrinsic time delay 
originating in the lens geometry. We discuss the feasibility of observationally distinguishing 
these two sources of time delay. We approach this issue through application of Fermat's 
principle, a useful tool for analysing gravitational lens problems |||§, which has recently 
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been shown to hold in the non-stationary space-times we consider . We conclude that 



measurements of time delays in gravitational lenses are not likely to serve as a method to 
detect or constrain a cosmological background of gravitational waves, because the wave- 
induced time delay is observationally indistinguishable from the intrinsic time delay of an 
alternative lens geometry. As a consequence, the cosmological applications of lens time de- 
lays, e.g., inferring Ho or galaxy masses, are not affected by gravity waves, regardless of 
the amplitude h. We note that, using quite different methods, the same conclusions were 
reached for general (scalar, vector, and tensor) metric perturbations by Frieman, Kaiser, 



and Turner. 12 



In Refs. Jill, the time delay induced by a gravitational wave upon a gravitational lens 
was evaluated through the Sachs- Wolfe formula [|IJ for the differential photon redshift in the 
presence of metric fluctuations, integrated along unperturbed photon paths, i.e., along the 
same trajectories the photons would have followed in the absence of the wave. As we will 
show, this method is not applicable in the case that the wave amplitude h is comparable to 
or larger than the angular separation 2rj that the images would have in the absence of the 
wave, and the expression for the time delay derived in Refs. QJ^j is valid only if h <C 77. In 
the opposite limit, h 3> rj, the effect of the wave is equivalent to a change in the alignment 
of the system so large that multiple images do not form (at least for non-singular lens 
potentials). Thus, the wave- induced time delays never exceed typical intrinsic delays, and 
cannot be used to constrain the amplitude of cosmological gravitational waves. Moreover, 
even in cases where the wave-induced delay is comparable to the typical intrinsic delay, we 
will show that an observer would attribute the entire delay to the intrinsic lens geometry. 
Thus, the wave-induced delay cannot be unearthed in practice or in principle. 

To address these issues, we explicitly take into account the spatial distortion of the photon 
trajectories induced by the gravitational waves, which is non-negligible even if h << i]. 
The wave-induced perturbation of the photon paths gives rise to extra contributions to the 
time delay, in addition to the differential redshift along the two trajectories. The extra 
contributions arise as a consequence of difference in path lengths and different gravitational 



potential traversed by each photon due to the asymmetry in their trajectories induced by 
the wave. When the dust settles, our result for the wave-induced time delay coincides with 
that of Refs. because these extra terms cancel each other, but only in the limits h <^7] 
and uLrj <C 1. Moreover, the spatial distortion of the photon paths is always very significant 
when it comes to the interpretation of lens observations: a gravitational wave distorts the 
apparent angular positions of the images relative to the deflector in just such a way that an 
observer would attribute the wave-induced time delay to an intrinsic time delay associated 
with the image-deflector misalignment he or she sees. Since the lens geometry is not known a 
priori, but reconstructed from observations, one could equally well adjust the measurements 
to a given lens geometry in the presence of gravitational waves, or to an alternative lens 
geometry and no waves at all. Thus, it appears observationally impossible to distinguish 
wave-induced time delays from intrinsic delays, and so to detect cosmological gravitational 
waves through time delay measurements in gravitationally lensed quasars. 



II. TIME DELAY IN A SIMPLE LENS CONFIGURATION 

To more clearly display the features discussed above, we first analyze a simple lens model: 
a Schwarzschild (point mass) lens in a highly symmetric configuration, and a gravitational 
wave propagating perpendicular to the lens axis. In the next section we generalize the 
results derived here to the case of an arbitrary thin lens and arbitrary polarization and 
wave vector of the gravitational wave. Consider a static, spherical body of mass M, located 
at the origin of coordinates, that deflects photons emitted by a point-like source located 
at (x = 0, y = 0, z — —L), and an observer on the extension of the source-deflector line at 
(x = 0,y = 0,z = +L) (see Fig. 1). Given the axial symmetry, the observer sees an Einstein 
ring image of the source, with angular radius rj = \jlGMjL. Here, G is Newton's constant, 
we take the speed of light c = 1, and we assume rj << 1. To simplify the discussion, we focus 
on those photons that travel along the y = plane, forming two images on opposite sides 
of the ring. In the absence of a gravitational wave, there is no time delay between the two 
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images, and they arrive with an angular separation A8 = 2r]. Now consider a gravitational 
wave of dimensionless amplitude h and frequency lu, with polarization (+), propagating 
along the x axis in the positive x direction. Sufficiently far from the deflector mass, the 
spacetime interval can be approximated by 

ds 2 = (l - ^— ^ dt 2 - (l + 'L-^£\ (fa 2 + dy 2 + dz 2 ) + h cos u{t - x)(dy 2 - dz 2 ) . (2.1) 

Along a photon path, ds 2 = 0. Thus, if the spatial photon trajectories x = x(z) were known, 
one could evaluate the time of travel by simple integration in z from — L to L, 

T m J L dz 

To the level of approximation we shall be working (we are interested in terms of order hrj 
in the time delay), t can be replaced in eq. ( |2.2|) by t = t e + (z + L), with t e the time at 
which the photons were emitted at (x — 0, z — —L). The first two terms in the integrand of 
eq. ( |2.2| ) are the geometric contribution to the time of travel, while the third and fourth are 
contributions from the gravitational potential of the wave and the deflector respectively. 

A. Integration along unperturbed paths 

Let us first evaluate the time delay along unperturbed photon trajectories. We approxi- 
mate each path by straight segments, deflected by angle a = 2r] at the deflector plane z = 
(see Fig. 1). The approximation by straight segments is convenient and appropriate in 
the case of thin gravitational lenses, where most of the deflection occurs in the immediate 
vicinity of the deflector plane 0. The light trajectories are then 

x 1>2 = ±i](z + L) z < 

xi, 2 = Tv(z -L) z>0 (2.3) 

where the subscript (1,2) distinguishes trajectories that pass along opposite sides of the 
deflector. Straightforward integration leads to the time of travel, T\^. Clearly, given the 



1 I ' dx\ 1 



2GM 



(2.2) 
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symmetry of the integration paths, the only contribution to the time delay comes from the 
different gravitational wave potential encountered by each trajectory, i.e., from the third 
term in the integrand of eq. ( |2.2| ) . The time delay is 

IlTl 

AT = T\ — T 2 = [sinix>(t e + 2L) + sinujt e — 2sinuj(t e + L) cos u Lrj] , (2-4) 

which, if uoLt] « 1, is approximated by 

AT « 4^ sin 2 (^) sin u{t e + L) . (2.5) 
to V 2 J 

This coincides with the result of Ref. |3J], evaluated by integration of the Sachs- Wolfe formula 
along the same unperturbed trajectories. 

The method used above to evaluate the time delay induced by the gravitational wave is 
questionable, even if h « rj, because the actual photon trajectories are perturbed by the 
wave, and are expected to be neither straight nor symmetric with respect to the lens axis. 
As a result, one path may have smaller impact parameter with respect to the deflector than 
the other, and hence be deflected by a larger angle. This asymmetry in the paths leads to 
differences in both the geometric and potential contributions to the time of travel as large 
as that evaluated above. We shall evaluate these extra contributions, after derivation of the 
lens equation through Fermat's principle, and find that they cancel to leading order only if 
h <C 77. Thus, eq. ( |2.5|) is only valid in this limit (moreover, in this limit, the wave-induced 
delay would be swamped by the intrinsic delay if the perfectly aligned symmetric lens were 
replaced by a misaligned lens such as that in Fig. 2). Furthermore, we will show that, even 
in this limit, the delay (|2.5|) would be attributed by the observer to intrinsic lens delay, that 
is, to an apparent misalignment between source, deflector, and observer. 

To prepare for this result, we first briefly review the derivation of the lens equation and 
time delay for a misaligned lens in the absence of gravity waves. 
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B. Fermat's principle for a Schwarzschild lens 



Fermat's principle provides a useful shortcut in many lensing problems: one can approxi- 
mate the photon trajectories by null zig-zag trial paths and then extremize the time of travel 
HII, instead of solving the geodesic equations. Consider the lens depicted in Fig. 2. The 
source is at an angle (3 with respect to the line that joins observer and deflector, and in this 
subsection we assume there is no gravitational wave present. Consider a path that, starting 
from the source at z — —L, moves along a straight line up to the z = plane, where it is 
deflected by an angle a, and then arrives at the observer at z = L, forming an angle 9 with 
respect to the line that joins observer and deflector. The angles a, (3 and 9 must satisfy 

a = 2(9 -(3). (2.6) 

Assuming (3,9 « 1, the time of travel along such a null path would be 

T w 2L + 9 2 L - 2/39L - AGM In 9 , (2.7) 

where we have neglected constant (^-independent) terms. The first three terms are of purely 
geometric origin, while the last originates in the gravitational potential of the deflector. The 
condition that T be an extremum (dT/d9 = for fixed (3) gives 

9*-(39- 2 -^ = 0, (2.8) 

which is usually referred to as the lens equation |§ . The solutions give the angular positions 
of the two images: 

e lt2 = ±(J3±JjP + 4tf) with 7? 2 = ^. (2.9) 

Different signs indicate that the images appear on opposite sides of the deflector. The 
resulting time delay is given by 

AT = T 1 -T 2 = (9\ - 9\)L - 2/?(0 a - 9 2 )L - 4GMln(0i/|0 2 |) • (2.10) 

In the limit of small misalignment angle, (3^7], the image angular positions are approxi- 
mately 
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6 1 , 2 tt±r ] + ~(3 . (2.11) 

and the time delay reduces to 

AT « -4(3t]L » -2/3A0 , (2.12) 

where A# = 61 — 62 is the image angular separation. 

For the arguments we will make below, it is useful to bear in mind how lens observations 
in such a simple system could be used to extract cosmological information. If the deflector is 
seen in addition to the two images, then the lens observables are 61, 62, and AT. The observer 
can then infer the misalignment angle from f3 = 9\ + 6 2 , and the lens parameter 7] from eq. 
Q2.9|). Using these observed and derived quantities, eq. (|2.10| ) can be used to determine L. 
(More generally, if the source and observer are not equidistant from the lens, the reasoning 
above determines a distance measure for the lens.) Comparison with the deflector redshift 
then yields an estimate of the Hubble parameter Hq. For deflectors more complex than point 
masses, the observables above must be supplemented by information about the deflector 
potential obtained, e.g., from measurement of the lens' velocity dispersion. 



C. A Schwarzschild lens with a gravitational wave 

Now we proceed with a similar technique, based on Fermat's principle, to evaluate the 
time delay induced by a gravitational wave in the symmetric Schwarzschild lens configu- 
ration, with source, deflector, and observer aligned as in Fig. 1. The validity of Fermat's 
principle in non-stationary spacetimes was recently discussed in the context of gravitational 
lensing problems [ 10 , . We evaluate the time of travel, integrating eq. fl2.2| ). Instead of 



integrating along straight lines, however, we take the two segments of each zig-zag trial path 
to be null geodesies of the gravitational wave metric, as they would have been in the absence 
of the lensing body. We work up to the order of approximation needed to study terms of 
order hrj in the time delay, and we also assume uLrj << 1. This is the most interesting 
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range, since the effect under consideration becomes largest when ujL pa 1. The geodesic 
equations in the metric ( |2.1| ), with M — 0, lead to 

dx 1 

— pa 7 - -hcosuj(t e + z + L) . (2.13) 
2 

Here, 7 is an arbitrary integration constant, the average slope of the trajectory, assumed to 
be small. Note that, as before, in the argument of the cosine in (|2.13 ), we have replaced t by 



t e + z + L; we have also dropped the dependence on x because it is unnecessary to include 
it to evaluate the time delay to order hi], in the limit ujLt] << 1. The third term in eq. 
fl2.2Q is the only one where the x-dependence inside the cosine needs to be included, and it 
is enough to do so at zero order. 

Now we choose the integration constants so that a trajectory that starts from x = 0, z = 
—L at t — t e , and deflected by an arbitrary angle in the z = plane, arrives at the observer 
at x = 0, z = L. One finds that 

dx I 

— = e hcosu(t e + z + L) if z < 

dz 2 

dx h I 

— = -e-\ [sin uj (t e + 2L) - sin ujt e ] h cosu(t e + z + L) if z>0. (2.14) 

dz 2ujL 2 

Here, e is an integration constant which parametrizes the family of trajectories that meet the 
focusing conditions at the required points. Each trajectory consists of two segments which 
are null geodesies of the gravitational wave metric, neglecting the deflector potential, which 
will be taken into account through Fermat's principle. According to the latter, the actual 
trajectories are those null paths that extremize the time of travel with respect to variations 
of the parameter e in the metric that includes both the gravitational wave as well as the 
deflector's potential. 

Before we proceed to extremize, however, we parametrize the trajectories in a different 
way, defining a parameter more relevant to observations. Instead of e, we use the angular 
position of the image, which we denote by 9, relative to the angular position of the deflector, 
at the time of arrival of the images at the observer. We again use equation ( j2.13j ) and fix 
the appropriate integration constants to determine the slope of the trajectory of a photon 
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that arrives from the deflector (i.e., the angular position of the deflector), which we denote 
by dx lcns /dz, 



dxi e 



dz 



= - — -[sm u(t e + 2L) — smu(t e + L)] — -hcosto(t e + 2L) 
-l 2uL 2 



(2.15) 



Then the angular position 9 of the image with respect to the apparent deflector position is 

h 



dx 
dz 



dx\ en 
=l dz 



■[smut e — sincj(t e + L)] . 



(2.16) 



=i ' 2ujV 

This relates 9 to the parameter e of eq. ( 2.14Q . Now the deflection imprinted by the lens 
upon the trajectory at z — 0, which we denote by a, can be written as 



a 



(l.i 

dz 



2=0" 



dx 
dz 



=0+ 



2(6 -fig 



(2.17) 



where we have defined 
h 
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[sinc<j(t e + 2L) + sinc<jt e — 2sincj(t e + L)\ 



h . 2 / toL\ . 

sir si 



T sin , 

ujL V 2 



\ smu(t e + L) . 



[2.1%) 



We have defined the quantity (3 g in such a way that the expression (|2.17|) for the deflection 
has the same form as eq. fl2.6|) -in that case, f3 measured the misalignment between deflector, 
source, and observer in the absence of a gravitational wave, as in Fig. 2. We will see in 
what follows that in all respects j3 g plays exactly the same effective role here. 



Next we evaluate the time of travel, integrating eq. ( |2.2| ) along the null trajectories 
( j2.14| ), parametrized in terms of 9, and find 

T^2L + 9 2 L-2(3 g 9L-AGM\n9. (2.19) 

As advertised, this has exactly the form of eq. ( |2.7|) , which gave the time of travel for a 
similar lens with no gravitational wave, but with lens and source misaligned by an angle 
(3, as in Fig. 2. Recall that in that case, the first three terms were of geometric origin. 
In the present case, only the first two terms are geometric (they come from integration of 
[l + (dx/dz) 2 /2] in ( |2.2| )). The third term, proportional to (3 g , is due to the wave gravitational 
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potential, and comes from integration of the third term in eq. ( |2.2|) . Finally, the last term 
is due to the deflector's gravitational potential. 

The equivalence of expressions ( |2.7|) and ( |2.19| ) leads to our main conclusion: the lens 



equation in the presence of a gravitational wave is, to the order of approximation considered, 
completely equivalent to that of a similar lens with a different alignment and no gravity wave. 



The effective misalignment angle (5 g is given by eq. ( |2.18|) in terms of the wave parameters. 
The analogy is exact only in the limit uLr] -C 1, but this is the interesting range in any case, 
because 77 1 and time delays are largest for uoL m 1. Since f3 g depends upon time, the 
analogy is only valid over periods of time much shorter than oj~ 1 . Again, since the effect is 
relevant only for waves of cosmological wavelengths, this time- variation of the time delay is 
observationally irrelevant. 

From eq. ( |2.19|) , the time delay between the two images is given by expression ([2.1 0|) 
with the substitution j3 — > (3 g , 

AT = T 1 — T 2 = (61 - 9l)L - 2p g {6x - 6 2 )L - 4GMln(0i/|0 2 |) . (2.20) 

So far, we have made no assumption about the relative amplitudes of the gravitational wave 
effect, (3 g ~ h, and the deflector Einstein ring angular radius 77 (Cf. eq. ( |2.9| ) ). However, if 
P g 3> f], the effect of the wave is equivalent to that of a system very much out of alignment. 
In this limit, the magnification of the second image goes to zero as {fj/Pg) 4 , and multiple 
image formation effectively does not take place. 

In the opposite limit, (3 g <C 77, Fermat's principle leads to the same result as eq. fl2.12|) , 
but with (3 replaced by the effective (3 g of eq. ( |2.18|) , 

AT w -A/3 g r]L = 4— sin 2 (—) sin uo{t e + L) . (2.21) 



u \ 2 

This result coincides with that of eq. ( p.5| ), obtained through integration along unperturbed 
paths, and is just Allen's result [Q. Note that the additional term originating in a path 
length difference, (Of — 6%)L, cancels the term due to the deflector's gravitational potential, 
AGM ln(#i/|#2|)- The wave-induced distortion of the photon paths is not negligible, however, 
when it comes to interpreting the result. 
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Indeed, suppose the observer of this lens has no knowledge of the possible existence 
of gravitational waves, and seeks to measure, e.g., the deflector mass M or the Hubble 
parameter Hq from her observations. The effect of the gravity wave upon the apparent 
angular positions of the images and the deflector would trick the observer into believing he or 
she sees an ordinary misaligned lens. Moreover, the observer's inference of the misalignment 
angle (3 g from the observed image angular positions, (3 g = 9\ + 9 2 , and the observed image 
time delay would all be in accord with this belief, and he or she will infer the correct values 
for M and Hq, even though taking no account of gravitational waves and instead assuming a 
homogeneous and isotropic spacetime (aside from the deflector). That is, while the gravity 
wave does cause a time delay, it covers its own tracks in a misalignment change, leaving 
no measurable trace of its presence, and can be safely and consistently ignored by the lens 
observer. Thus it appears impossible to use time delay measurements to detect cosmological 
gravitational waves even in principle. 

For completeness, we emphasize that this conclusion holds even if /3 > r), but that it is 
only in the limit (3 g -C 77 that the time delay agrees with eq. ( |2.5| ). 

III. TIME DELAY IN A THIN LENS WITH GRAVITATIONAL WAVES 

In this section we show how the conclusions reached above can be generalized to the 
case of an arbitrary thin gravitational lens with gravitational waves of arbitrary polarization 
and direction of propagation. First we briefly review the features of a general lens when 
no gravitational waves are present. We assume, as usually applies for cases of astrophysical 
interest, a thin, stationary gravitational lens, such that the weak field approximation is valid 
||. Consider a lens geometry as in Fig. 2, only now we do not assume that the photon paths 
lie in the plane that contains source, deflector, and observer: 9 and (3 are two-component 
angular vectors, that give angular positions at the observer's location, a is the deflection, 
and £ = L9 determines the impact parameter in the deflector plane. The condition that 
the photons from the source reach the observer implies the following relation, the vectorial 
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generalization of eq. 

a = 2 (9-0), (3.1) 

and the time of travel is given by: 

T^2L+\e\ 2 L-2(3-6L-tlj(0, (3.2) 

which is the generalization of eq. (|2.7|) . The last term originates in the deflector gravitational 
potential, and is given, for a thin lens, by 

iKQ =4Gj d 2 e £(f) In (^j^) , (3-3) 

where £ is the mass density projected on the lens plane. Note that this term depends only 
upon the impact parameter £, reflecting the fact that in the thin lens approximation the 
effect of the gravitational potential of the lens is dominated by that part of the trajectory 
closest to the deflector. Here £o is an arbitrary length scale. Following Fermat's principle, 
we extremize the time of travel with respect to 9 and arrive at the lens equation: 

*-?-±% = 0. (3.4) 

Notice that dip/ '89 = Ldip/d^. The solutions to this equation give the angular positions of 
the images. 

Now we show the equivalence between the effect of a gravitational wave and an effective 
lens misalignment. Consider a lens geometry with deflector, source, and observer aligned at 
z = 0, —L, L respectively, as in Fig. 1, along the z-axis. Let U be the gravitational potential 
of the deflector. Consider a gravitational wave propagating at an angle •d with respect to 
the lens axis. We take the (x, z) plane as that containing the lens axis and the direction of 
propagation of the gravitational wave. The metric perturbation caused by the wave can be 
expressed as: 

^ —cos 2 dh+ — cos^/iv sin^cosi? h ^ 



hij 



— cos$ h x h + smd h x 

sin cost? h+ sin^ h x — sin 2 $ h 



cos(ut — k ■ x) (3.5) 
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with propagation vector k = u;(sin$, 0, cos^), and h + and h x the amplitudes of the two 
wave polarizations. The total metric is then given by: 

ds 2 = (1 + 2U)dt 2 - (1 - 2U){dx 2 + dy 2 + dz 2 ) + h lj dx i dx j . (3.6) 



Along a null path, the time of travel is given by 



7 % [ L dz 



1 ( dx\ 2 1 / dy\ 2 1 dx % dxi 

i+ 2u + 2 y + 2^ii7- 2c/ 



(3.7) 



Notice that / 2U dz is the same as what we had previously defined by ip in eq. ( |3.2|) . 

We now define a family of null trial paths along which we will integrate eq. ( |3.7|) . Each 
path is built out of two segments deflected by an angle a at the deflector plane. Instead of 
taking straight trajectories, we let each segment be a solution of the geodesic equations in 
the presence of the gravitational wave, neglecting the potential U of the deflector, since its 
effects are later taken into account through Fermat's principle. The condition that a photon 
from the source at (x, y = 0, z = —L) reaches the observer at (x, y — 0, z — L) defines a 
one-parameter family of trajectories parametrized by an arbitrary vector e = (e x , e y ): 

— = e x — sini?(l — cos'i?) cosuj(t e + L + z(l — cos^)) 

dz 2 

dy 

— = e v + h x smi!) cosoo(t e + L + z(l — cost?)) if z<0 , 
dz 

—= —e x + — y- sini?[sina;(£ e + L(2 — cos^)) — smu(t e + L cost?)] 
dz 2oo Li 

h 

— — sini?(l — cos'i?) coscj(t e + L + z(l — cos^)) 



dy h x sin^ 



— e„ - 



-[sm u(t e + L{2 — cos$)) — smu(t e + L cosi?)] 



dz " y cjL(1-cosi?) 

+h x sin$coscu()!; e + L + z(l - cos$)) if z>0. (3.8) 

The wave also affects the apparent position of the deflector; as before, we change variables 
from <f to the relative angular position between the image and the deflector at the observer's 
position, which we denote by 9. We find the relation (assuming ujLO « 1): 

n f h + sm$ —h x sin^ \ 

= {ex,€y) + T — , —rr, [smu(t e + LcosV) - sm u(t e + L)\ . (3.9) 



2ooL ' uL{l — cos?9) 
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The deflection a imprinted upon the trajectory at the lens plane z = can be written as 



a = 2(6- fa) 



(3.10) 



if we define j3 g as 



Po = 



h. 



-sin??, 



sin-;? 



2uL (1 -costf), 
x[sinu;(t e + L(2 — cost?)) + sinu;(£ e + Lcos-#) — 2sinu;(t e + L)] 



' h+ . —h x 2 sin 
uL ' tuL (1 — cost?) 



sm 



(1 — cost?) 



sina;(t e + L) 



(3.11) 



and (|3.11|) reduce to our previous eqs. (|2.16| ) and (|2.18|) in the case 



Of course, eqs. 
= n/2,h x = 0. 

Now we are ready to find the time of travel, integrating eq. (|3.7|) . One important thing 
to note is that, for a thin lens, integration of the last term in ( |3.7|) , the contribution of the 
deflector gravitational potential U, gives — where £ = x(z = 0) is the impact parameter 
of the trajectory. And, if uL6 « 1, the relation £ = L9 still holds, as in the absence of a 
gravitational wave. At the end, we find for the total time of travel: 



T = 2L+ \9\ 2 L - 2 ■ 9 L - ip(L6) 



(3.12) 



Since this has exactly the same functional dependence on 9 as in eq. (|3.2|) , we confirm the 
equivalence between an aligned lens in the presence of a gravitational wave and a lens with 
an effective lack of alignment given in terms of the wave parameters by (3 g of eq. ( |3.11|) . 

In the special case of an axially symmetric lens and in the limit \(3\ « \9\, the solutions 
#*i 2 to the lens equation lie in the plane that contains the lens axis and the direction of j3 g . 
That plane forms an angle with the plane that contains the gravitational wave propagation 
vector (which we took as the (x, z) plane) such that 

(Pg)y 2/i x 



tan ( 



{Pg) x h + (l-COSti) 



(3.13) 



The solutions are then of the form 



9 1;2 ~ ±ff + aP 
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(3.14) 



with 



77 = t/(cos0, sin 1 



V 



1 dijj 
2L~d6 



e=rj 



(3.15) 



a solution to the unperturbed lens equation, and a a coefficient that depends upon the lens 
model: 



-1 



1 d 2 ip 



8=r) 



(3.16) 



2L 89 

For instance, a = 1/2 for a Schwarzschild lens, where ip = 4GMln6, and a = 1 for a singular 
isothermal sphere, where ip = 47rE^#, with E„ the velocity dispersion. For a non-singular 
lens, there must be an odd number of images; in that case, eq. (|3.14|) refers, say, to the 
outer two images (the central third image is usually de-magnified). 
The time of travel for these solutions can be expanded as 

dip 



T 1)2 = 2L + (\ff\ 2 ± 2a(3 ■ ff)L T^P-fjL- ip(±Lff) - a(3 



09 



(3.17) 



Using the unperturbed lens equation ( p,15[ ) we see that the two contributions to the time of 
travel proportional to a/3, one of geometric origin and the other due to the deflector's grav- 
itational potential, cancel each other. Besides, ip{Lff) = ip(-Lff) for an axially symmetric 
lens. Finally, the time delay between two images for a thin, axially symmetric lens, in the 
limit \j3\ « \ff\ and with loLt] « 1 is 

2 sin 6 



4r? 

AT w -A/3-ffL = — sin$ 
to 



h + cos0 + h y 



costf) 



sm 



— cosfi) 



sino;(t e + L) 



(3.18) 



where <fi is the angle between the plane containing the photon trajectories and the plane 
that contains the gravitational wave propagation vector, as given by ( |3.13| ). In the limit 
(3 g « 77, expression ( |3.18| ) agrees with that of Ref. M, evaluated through integration of the 



Sachs- Wolfe formula along unperturbed photon paths. 



17 



IV. CONCLUSIONS 



We have shown that the lens equation for a thin, axially aligned gravitational lens con- 
figuration in the presence of a very long wavelength gravitational wave is equivalent to that 
of a similar lens with the source out of alignment and no gravitational wave. An observer 
who measures time delays, angular positions, or any other observables such as relative mag- 
nifications and redshifts, and uses them to reconstruct the lens configuration, cannot tell 
the two situations apart. Thus, an observer ignorant of gravitational waves would naturally 
and 'correctly' interpret the observations as a simple non-aligned lens. This conclusion is 
valid if ujLti « 1, which is the interesting range since the induced time delays are largest 
when uoL Pd 1. We performed our calculations around an aligned lens configuration and with 
source and observer equidistant from the deflector, but it is clear that the conclusion holds 
in more general cases: the effect of a long wavelength gravitational wave upon a lens with 
a given geometry is equivalent, from the observer's viewpoint, to an effective change of lens 
geometry. Consequently, measuring time delays in gravitational lenses does not provide a 
method for probing a cosmological background of gravitational waves. 

Formally, Allen's result for the wave-induced time delay is correct in the small amplitude 
limit: for (3 g << rj, eq. (|3.18 ) for the time delay induced by a gravitational wave upon a 
thin, axially symmetric lens agrees with that of Ref. ||. In the opposite limit, (3 g ^> 77, the 
effect of the gravitational wave is equivalent to a change in the alignment between source, 
deflector, and observer by an amount that exceeds the typical deflection angle the deflector 
can imprint, precluding the formation of multiple images in the case of an aligned lens. In 
this case, multiple images can only form if there is a compensating geometric misalignment 
between source and deflector, and the geometric delay will partially cancel the lens-induced 
delay. Thus, even ii h » r], the total time delay does not exceed the typical intrinsic lens 
time delay of order Lrf. Moreover, in either limit, the measured time delay is just what 
the observer would expect in the complete absence of gravitational waves, based on her 
measurements of the lens observables. 
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The detection by the COBE satellite [Q] of a quadrupole anisotropy in the cosmic mi- 
crowave background places a bound ft<6x 10~ 6 on the amplitude of cosmological gravita- 
tional waves. If a large fraction of the anisotropy detected by COBE is due to gravitational 
waves, a possibility that can be accommodated by many inflationary cosmological models 
then the wave-induced time delays between multiple images of quasars are comparable 
to typical intrinsic lens time delays, with rj « 10~ 5 . One could have hoped that careful 
lens modelling could allow one, at least in principle if not in practice, to separate the wave- 
induced from the intrinsic time delay, and thus reveal the presence of the gravitational waves. 
Our work indicates that this is not a possibility. 
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X 




Figure 1 : The geometry of an aligned lens. The deflector is a point mass M at the origin 
of coordinates. Source and observer lie along the z-axis, equidistant from the deflector. 
The observer sees two images of the same source (actually an Einstein ring) with angular 
separation 2i]. The trajectories are approximated by straight segments. 
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Figure 2 : The geometry of a non-aligned lens. The source forms an angle (3 with respect 

— * — * 

to the line that joins observer and deflector. 9 1 and 9 2 are the angular positions of the 
images, £ is the impact parameter, and a the deflection angle. 
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